Coherence Resonance in an Autapse Neuron Model with Time Delay 
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We study the noise activated dynamics of a model autapse neuron system that consists of a 
subcritical Hopf oscillator with a time delayed nonlinear feedback. The coherence of the noise 
driven pulses of the neuron exhibits a novel double peaked structure as a function of the noise 
amplitude. The two peaks correspond to separate optimal noise levels for excitation of single spikes 
and multiple spikes (bursts) respectively. The relative magnitudes of these peaks are found to be a 
sensitive function of time delay. The physical significance of our results and its practical implications 
in various real life systems are discussed. 
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The constructive role of noise in the dynamics of com- 
plex systems is a subject of much current interest and ac- 
tivity. Important manifestations of such a behaviour are 
seen in basic phenomena like stochastic resonance (SR) 
[l|, 13, I31j coherence resonance(CR) 4] or noise-induced 
synchronization of dynamical systems |5|, |6J . Coherence 
Resonance (CR) which refers to the resonant response 
of a dynamical system to pure noise, is closely related 
to the phenomenon of stochastic resonance and is some- 
times also known as Autonomous Stochastic Resonance 
(ASR) 7\. The effect, first noticed by Sigeti and Hors- 
thcmke 8] in a general system at a saddle-node bifurca- 
tion, implies that a characteristic correlation time of the 
noise-excited oscillations has a maximum for a certain 
noise amplitude. This has been clearly demonstrated for 
the classic FitzHugh-Nagumo neuron model 4j and shown 
to have a deep connection to the excitable nature of the 
system. CR can have important consequences for neuro- 
physiology or other complex systems where a significant 
degree of order can arise through interaction with a noisy 
environment. Past studies of CR have been mainly con- 
fined to simple systems whose noise-induced nonlinear 
outputs consist of impulsive excitations of a single kind 
e.g. spikes which have two characteristic time scales - a 
fast rise time and a longer decay time. This phenomenon 
has been found not only in various lab experiments, such 
as electronic circuits Q, laser systems 10, IJ, electro- 
chemistry [l2|, or BZ reactions [13 but also in natural 
systems such as ice ages in climatology [lj| or dynamos 
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As is well known, excitable systems especially in 
neuroscience, can also have a more complex response 
in terms of various time scales, such as short time 
spikes and bursts (multi-spiking) with different tempo- 
ral signatures fla \11\ . The nature of CR in the presence 
of different kinds of excitations (e.g. spikes and bursts) is 
not known and is one of the important objectives of the 
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present study. In order to explore this question, we study 
the noise activated dynamics of a novel model neuron 
system that was recently presented in [ig. The model 
consists of a subcritical Hopf oscillator with a time de- 
layed nonlinear feedback. It contains the necessary in- 
gredients to reproduce most of the basic features of ex- 
citability that are displayed by standard models such as 
the Hodgkin-Huxley system |l9l | and in addition provides 
a convenient means of studying the effect of time delay 
on neural dynamical behavior. The basic mathematical 
form of the model is, 

Z{t) = [i{uj + b\z{t)\^) + |z(t)|2 - |z(t)|4] z{t)-kz\t--T) 

(1) 

where z = x + iy is a complex amplitude and the fre- 
quency of the oscillations is determined by uj and b\z\'^. 
The parameter b which is called the shear parameter, de- 
termines how the frequency depends on the amplitude 
of the oscillations, k is the magnitude of the feedback 
strength and r is the time delay parameter. In the ab- 
sence of the feedback term the oscillator is poised at 
the subcritical bifurcation point. The nonlinear feed- 
back term, which provides the basic excitable behavior, 
is modeled to mimic the structure of the recently discov- 
ered autapse neurons which have axons terminating back 
on their own dendrites ^20l |21| . The feedback is time de- 
layed to account for finite propagation times of signals. 
The basic dynamical behaviour of the system is best cap- 
tured in a two parameter bifurcation diagram (in k and 
T space) which was obtained in 18] and is reproduced 
here as Fig. ^ The diagram delineates the different bi- 
furcation branches as well as the regions of stable limit 
cycle, stable fixed point and bistability. When subjected 
to an external periodic signal and noise the system dis- 
plays both spike trains as well as multi-spiking (bursty) 
behaviour [l^. For our present CR studies, we examine 
the temporal response of the neuron in the presence of 
only an external noisy stimulus, namely /(t) = \/2D^{t), 
where ^(i) is zero mean Gaussian white noise with inten- 
sity D. We confine ourselves to values of k that are above 
the critical value of kc — 0.42506 and to r values below 
the bi-stable region and choose different noise strengths 
D. 
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FIG. 1: Stability diagram in the parameter space of k and 
r. Note the various bifurcation boundaries demarcating the 
different stability regions. (SNLC; Saddle-Node on Limit Cy- 
cle bifurcation, SN: Saddle-Node bifurcation, SSL: Saddle- 
Separatrix Loop bifurcation) 
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FIG. 2: The CR measure 7? and mean interspike interval {T) 
versus the noise intensity D. The solid curve is for r = and 
the dashed curve is for t = 0.30. The feedback strength is 
k = 0.426. The point is close to SNLC bifurcation branch. 



For each noise intensity D we execute rather long sim- 
ulations (100 datasets consisting of 200,000 time steps 
each with St — 0.05) and collect a large number of inter- 
spike intervals (ISI) T. Using this data we determine the 
standard parameter for coherence resonance i?, which is 
given by the ratio of the mean of the int ersp ike intervals 
((T)) and its standard deviation (ctt )|l[l7[l2ll2^. 



R={T)/aT 



(2) 



If the data consists of a completely random (Poissonian) 
set of spikes then R would have a value of unity, whereas 
a strictly periodic spiking train (maximum order) would 
make R — oo. In general R > I indicates a coherence, 
whereas values of R that are significantly lower than 
unity indicate lack of coherence i.e. random dynamics. 
Our results for the present neuron model are shown in 
Fig. 2 where the upper block shows the CR and the 
lower block depicts the average interspike interval (ISI) 
as functions of the noise intensity D. The solid curves 
are for r = and the dashed curves correspond to finite 
time delay (in this case r = 0.3). We observe that unlike 
the standard CR results of a unimodal response curve 3 
our present system produces a more complex response 
consisting of one well-expressed peak and, additionally, 
a broad shoulder in another region of noise intensities. 
The two-peaked structure is preserved in the presence 
of time delay with however one very substantial change 
- the relative amplitudes of the peaks are reversed. 
Without time delay there is a broadband peak for small 
noise intensities D around 0.003 and a shoulder with a 
minor peak at about D = 0.3, whereas in the presence of 
time delay (r = 0.3), we get a broad shoulder for small 
noise and a rather narrowband peak at D around 0.3. 
Note that in both cases the average ISI is decreasing 



with increasing D. Hence, the activity in the low noise 
region can be related to spiking (isolated spikes have a 
long interspike interval), whereas that in the strong noise 
region can be associated with bursting activity (several 
spikes within a burst leading to shorter ISIs). The 
bi-modal structure of the CR is a novel result which, to 
the best of our knowledge, has not been observed before 
and suggests that it can effectively distinguish between 
spikes and bursts and also provide a measure of their 
relative predominances. The presence of both spikes and 
bursts can be understood physically from the bifurcation 
diagram shown in Fig. 1. For r = and k — 0.426 the 
system is located on the bottom horizontal line in the 
region of a stable fixed point. The introduction of noise 
pushes the system randomly into the stable limit cycle 
region from which it returns into the quiescent region 
along an extended trajectory on the unstable manifold. 
For low values of noise intensity the system barely makes 
it into the limit cycle regime and the resultant outcome 
is a single spike. At larger values of noise the system 
is driven deeper into the oscillatory region and spends 
more time executing several limit cycle circuits before 
returning to the quiescent state in a spiked manner. As 
a result the temporal trace now consists of an onset 
spike followed by several oscillations and a terminating 
spike - all of which taken together constitutes a burst. 
Such a phenomenon is well known in the literature as 
the circle/circle or parabolic bursting [Ig. Spikes and 
bursts are optimally excited at two different values of 
the noise intensity as demonstrated by the CR curve. In 
the absence of time delay spikes appear to predominate 
as seen by the larger peak in the low noise regime and 
bursts have a relatively lower count of occurrences. 

With the introduction of time delay the relative 
predominance of bursting appears to grow as a function 
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FIG. 3: Probability distribution function of ISIs for k = 0.426, 
D = 0.03 and different values of r (= 0.00, 0.10, 0.15 and 0.30) 
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FIG. 4: Typical time series for k = 0.426 and D = 0.03 for 
(a) r = and (6) r = 0.3. Spikes dominate in (a) whereas 
predominantly bursting events are seen in (b) 



of time delay. This can be seen from the dashed curves 
in the CR curve of Fig. 2 where there is a relative rise in 
the second peak at higher noise values. To get a better 
measure of this increase in bursting activity, we have 
also looked at the normalized probability distributions 
of the ISI for various values of t at a fixed value of 
D = 0.03 (Fig. 3). The value of Z? = 0.03 is chosen 
because it marks a region in parameter space where for 
T — spikes and bursts are nearly equally abundant. 
This is clearly seen in the solid curve of Fig. 3, which 
has a two-hump distribution with the first peak (low 
ISI) corresponding to bursting and the second one of 
nearly equal height corresponding to spiking (large ISI). 
Note that the second peak also has a long tail which 
is characteristic of spiking. With increasing values of 
T, the second peak gradually diminishes (the various 
dashed curves) and the long tail gets suppressed. On 
the other hand, the first peak becomes substantially 
stronger and for t = 0.3, bursting behaviour clearly 
dominates. These results are also visually discernible in 
snapshots of segments of time series data corresponding 
to the r values of and 0.3 respectively in Fig. 4a, b. 

To understand the dynamical reason for this transition 
from spiking to bursting as a function of r, we have 
examined the effect of time delay on the characteristic 
properties of the limit cycle itself. We find that the 
primary influence is on the basic frequency of the limit 
cycle. A typical example is shown in Fig. 5 where the 
limit cycle frequencies for different values of the feedback 
strength k are shown for t = (solid circles) and t = 0.3 
(solid line). We see that the frequencies are always 
higher at finite values of t. The physical consequence of 
this change in frequency is that the system now executes 
a larger number of cyclic orbits before returning to 
the quiescent state - which is characteristic of bursty 
behaviour. The consequent preponderance of small ISIs 
as compared to the no delay case is directly reflected 



in the probability distribution function curves shown in 
Fig. 3. The sensitivity of the excitability property of the 
system to time delay is quite remarkable and is the other 
major interesting result of our present investigation. 
It suggests that time delay can provide an effective 
mechanism for steering the system towards either a 
predominance of spiking or of bursting behaviour. 
Since the nature of excitability is at the heart of many 
important neuronal functions such as communication, 
computational properties and information processing 
our finding can have important practical consequences 
for the collective dynamics of such interacting neurons. 
As an example the process of noise coupled synchro- 
nization between trains of neuronal signals can be 
significantly influenced by changes in the intrinsic time 
delay parameter of each neuron and create interesting 
consequences for inter-neuron communication. Likewise 
the computational properties of a neuron can change 
due to its transition from a spiking to a bursting state 
brought about by the presence of time delay. 

To conclude, in this paper we report two interesting 
results from the dynamical study of a noise driven 
model neuron that consists of a subcritical Hopf os- 
cillator with a time delayed nonlinear feedback. We 
find that the coherence resonance parameter of such 
a system possesses a bi-modal structure with the two 
peaks corresponding to optimal noise levels for two 
different kinds of excitations - namely single spikes and 
multi-spike (bursty) structures. Such a discriminatory 
ability enhances and to some extent generalizes the 
utility of the CR measure as a statistical tool for the 
analysis of neuronal data. Our other finding is that 
the presence of time delay in the nonlinear feedback 
can significantly influence the excitability properties 
of the system and bring about a controlled change in 
state from a predominantly spiking behaviour to largely 




FIG. 5: The frequencies of the hmit cycles for t — (sohd 
circles) and for r = 0.3 (solid line). Note the increase in 
frequencies for a finite value of time delay. 



bursting behaviour. This can have important practical 
consequences for various neuronal mechanisms (e.g. 
communication, information processing, computation 
etc.) where action potentials come into play. 



All the numerical simulations in the present paper have 
been carried out with the help of the software packages 
XPPAUT m and DDE-BIFTOOL j^. XPPAUT also 
provides an interface to the continuation software pack- 
age AUTO 26\. The authors thank C.S. Zhou for valu- 
able discussions. JK acknowledges the support via his 
Humboldt-CSIR research award. 
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